32nd International Cosmic Ray Conference, Beijing 201 1 




o 

(N 

> 

O 

On 
X 

6 

in 

o3 



> 
On 
O 
00 



X 



TeV gamma-ray survey of the northern sky using the ARGO- YBJ experiment 
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Abstract: The ARGO-YBI experiment is an extensive air shower array with full coverage RPC detectors located at 
Yangbajing (4300 m asl, Tibet, China). It is operated with high duty cycle (> 86%) and a large field of view (~ 2sr). 
It continuously monitors the entire overhead sky at 7-ray energies above 0. 1 TeV. In the talk, we will present the result 
of the northern sky survey (between declinations of -10° and 70°) from an analysis of 4 years of the ARGO-YBJ data 
(between July 2006 and February 2011). There are four known TeV sources observed with significance greater than 5 
S.D.. The significance from Crab Nebula is more than 16 S.D.. 90% confidence level upper limits to the flux from all 
directions in the sky are also presented, which vary from 0.09 to 0.53 Crab unit for Crab-like point sources. 
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1 Introduction 

In the last decade, great advances have been made in 
ground-based 7-ray astronomy and more than 100 very 
high energy (VHE) 7-ray sources have been observed. A 
new window using VHE 7-rays is progressively established 
to probe the non-thermal universe and the extreme physics 
processes in the astrophysical sources. The VHE 7-rays are 
the emissions of the relativistic particles, which are accel- 
erated by the astrophysical shocks that are widely believed 
to exist in the pulsar wind nebulae (PWN), shell type su- 
pernova remnants (SNR), active galactic nuclei (AGN) and 
so on. These shocks may accelerate protons or electrons. 
Relativistic electrons can scatter the low energy photons to 
VHE band, i.e. Inverse Compton (IC) emission, while rel- 
ativistic protons would lead to hadronic cascades and VHE 
7-rays are generated by the decay of secondary ir° mesons. 
Hence, VHE 7-ray observations are also important in un- 
derstanding the origin and acceleration of cosmic rays. 

Up to now, the number of VHE 7-ray sources is 107, which 
includes 61 Galactic sources and 46 extragalactic sources 
[1]. Most of the identified Galactic sources belong to PWN, 
SNR, and binary systems, however, about one third of them 
are still unidentified. The extragalactic sources are mainly 
composed of blazars, including BL Lac type objects and 
flat spectrum radio quasars. 4 radio galaxies and 2 star- 
burst galaxies are also observed. Due to the absorption by 
the extragalactic background light (EBL), which causes a 
substantial reduction of the flux, VHE 7-ray observations 
are limited to nearby sources. The most distant source to 
date is 3C 279 with z=0.536 [2]. 



Recent advances in VHE 7-ray observation are mainly at- 
tributed to the successful operation of image atmospheric 
Cherenkov telescopes (IACTs), such as H.E.S.S., MAGIC, 
VERITAS and CANGAROO, which made most of the dis- 
coveries when searching for counterparts to sources dis- 
covered at lower energies. To achieve an overall view of 
the universe in VHE 7-ray band, an unbiased sky survey 
is needed just like what has been done by Fermi and its 
predecessor EGRET at GeV 7-ray band, which detected 
1451 and 271 objects, including 630 and 170 unidentified, 
respectively [3, 4]. The H.E.S.S. collaboration has made 
great progress in surveying the inner part of the Galactic 
plane and discovered 14 7-ray sources [5], however, due to 
the limitation of small field of view (FOV) and low duty 
cycle, for IACTs is difficult to perform a comprehensive 
sky survey. As a result, although with sensitivity lower 
than that of I ACT, an extensive air shower (EAS) array, 
such as the Tibet AS7, Milagro and ARGO-YBJ experi- 
ments, is the only choice to perform a continuous sky sur- 
vey at VHE band. To date, several surveys of the VHE 
sky have been performed by AIROBICC in 2002 [6], Mi- 
lagro in 2004 [7] and Tibet AS7 in 2005 [8]. The latter 
two surveys resulted in the successful observation of 7-ray 
emission from Crab Nebula and Mrk421. The best upper 
limits at energies above 1 TeV are around 275^600 mcrab 
achieved by the Milagro experiment. In 2007, Milagro up- 
dated its survey of the Galactic plane and 3 new extended 
sources were discovered [9]. The ARGO-YBJ detector has 
a better sensitivity than any previous EAS array at 1 TeV 
and is continuously monitoring the northern sky with de- 
clinations from -10° to 70°. This paper presents the TeV 
7-ray survey of the northern sky using the 4-year data of 
ARGO-YBJ. 
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2 The ARGO-YBJ experiment 



The ARGO-YBJ experiment, located in Tibet, China, at 
an altitude of 4300 m a.s.l., is the result of a collaboration 
among Chinese and Italian institutions and is designed for 
VHE 7-ray astronomy and cosmic ray observations. The 
detector consists of a single layer of Resistive Plate Cham- 
bers (RPCs). 130 clusters (a cluster is composed of 12 
RPCs and each RPC is composed of 10 logical pads) are 
installed to form a carpet of about 5600 m 2 with an active 
area of ^93%. This central carpet is surrounded by 23 ad- 
ditional clusters ("guard ring") to improve the reconstruc- 
tion of the shower core location. The total area of the array 
is 110 m x 100 m. More details about the detector and 
the RPC performance can be found in [10, 11]. The high 
granularity of the apparatus permits a detailed space-time 
reconstruction of the shower profile and therefore of the in- 
cident direction of the primary particle. The arrival times 
of the particles are measured by Time to Digital Converters 
(TDCs) with a resolution of 1.8 ns. This results in an angu- 
lar resolution of 1.6 degrees for showers with N pa d — 20 
and better for increasing N pa d, tending to 0.3 degrees for 
showers with N pa d > 1000 [12]. In order to calibrate the 
18,360 TDC channels, an off-line method [13] has been 
developed using cosmic ray showers. The calibration pre- 
cision is 0.4 ns and the procedure is applied every month 
[14]. 

The central 130 clusters began taking data in July 2006, 
and the "guard ring" was merged into the DAQ stream in 
November 2007. The trigger rate is 3.5 kHz with a dead 
time of 4% and the average duty cycle is higher than 86%. 
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Figure 1 : The thick solid line is the distribution of signifi- 
cances from all directions on the northern sky map with de- 
clinations from -10° to 70° . The thin solid line represents 
the best Gaussian fit to the distribution of significances 
excluding those cells that have a distance to the Crab, 
Mrk421, MGRO J 1908+06 or MGRO J2031+41 shorter 
than 4 degree. 
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Source 


83.75 


22.15 


17.9 


Crab Nebula 


166.25 


38.25 


13.2 


Mrk421 


286.85 


6.55 


5.8 


MGRO J 1908+06 


307.85 


41.95 


6.3 


MGRO J203 1+41 








TeVJ2032+4130 



Table 1 : Location of all regions with an excess greater than 
5 S.D. 



3 Data analysis 

The ARGO-YBJ data used in this analysis were collected 
from 2006 July to 2011 February. The total effective ob- 
servation time is 1265 days. To achieve a better angular 
resolution, the event selections for 6 groups with differ- 
ent Np a d used in [15] are applied here and only events 
with zenith angle less than 50 degrees are used. The to- 
tal number of events after being filtered used in this work 
is 1.7 xlO 11 . No 7/proton discrimination is made in this 
analysis. In order to obtain a sky map, the celestial coor- 
dinates (right ascension and declination) are divided into a 
grid of 0.1° x 0.1° bins and filled with detected events ac- 
cording to their reconstructed directions. In order to extract 
an excess of 7-rays from each bin, the so-called "direct in- 
tegral method" [16] is adopted to estimate the number of 
cosmic ray background events in the bin. To remove the 
affection of large scale anisotropy, a correction has been 
applied which can be found in [15]. To take into account 
the PSF of the ARGO-YBJ detector, the events in a circu- 
lar area centered on the bin with an angular radius of ^70 
are summed together. The Li-Ma formula [17] is used to 
estimate the significance. 



4 Results 

The significance distribution from all directions with de- 
clinations from -10° to 70° is shown in Figure 1. There 
are many excess points with significance greater than 5 
s.d., which can be attributed to the four known VHE 7-ray 
sources, i.e. Crab, Mrk421, MGRO J 1908+06 and MGRO 
J2031+41 (or TeV J2032+4130). The distribution is con- 
sistent with the expectation from random background fluc- 
tuations after removing the regions within 4 degree from 
the four known sources. A Gaussian fit to this distribution 
has a mean of 6.5 x 10~ 3 and a ct=0.99. 

Figure 2 shows the significance map of the northern sky for 
the ARGO-YBJ data set. There are four sources clearly vis- 
ible in the map. Table 1 lists the location of all regions with 
significance greater than 5 S.D.. For each independent re- 
gion, only the bin with the largest significance is presented 
in Table 1 . The position is consistent with the ARGO-YBJ 
pointing error 0.2°. 

Recently, the Milagro experiment observed 14 of the 34 se- 
lected Fermi sources at a pre-trial significance of 3 S.D. or 
more at the representative energy of 35 TeV [18]. A simi- 
lar result is also obtained by the AS7 experiment [19]. As 
an experiment with lower energy threshold than Milagro 
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Table 2: ARGO-YBJ result for sources observed by Milagro with significance greater 3 S.D.. 



Name 


Type 


R.A. 


Dec. 


Milagro 


ARGO-YBJ 


AS7 


lev 


(OFGL) 




(deg) 


(deg) 


(S.D.) 


(S.D.) 


(S.D.) 


Assoc. 


J0534.6+2201 


PSR 


83.65 


22.02 


17.2 


17.8 


6.9 


Crab 


J0617. 4+2234 


SNR 


94.36 


22.57 


3.0 


0.7 


0.2 


IC443 


J0631. 8+1034 


PSR 


97.95 


10.57 


3.7 


0.7 


0.3 




J0634.0+1745 


PSR 


98.50 


17.76 


3.5 


1.8 


2.2 


Geminga, MGRO C3 


J 1 844. 1 -0335 




281.04 


-3.59 


4.3 


2.1 






J lyUU.U+UjJO 




ZoO.Ul 


j.yj 


i.b 


1 .0 


1 .U 




J1907. 5+0602 


PSR 


286.89 


6.03 


7.4 


5.5 


2.4 


MGRO J 1908+06 


J1923.0+1411 


SNR 


290.77 


14.19 


3.4 


2.0 


-0.3 


H.E.S.S. J1923+141 


J1954.4+2838 


SNR 


298.61 


28.65 


4.3 


0.7 


0.6 




J1958. 1+2848 


PSR 


299.53 


28.80 


4.0 


1.7 


0.1 




J2020. 8+3649 


PSR 


305.22 


36.83 


12.4 


0.9 


2.2 


MGRO J20 19+37 


J202 1.5+4026 


PSR 


305.40 


40.44 


4.2 


3.9 


2.2 




J2032.2+4122 


PSR 


308.06 


41.38 


7.6 


5.8 


2.4 


MGRO J203 1+41 


J2229.0+6114 


PSR 


337.26 


61.24 


6.6 


-0.8 




MGRO C4 


LAT PSR J2238+59c 


PSR 


339.56 


59.08 


4.7 


1.2 


2.5 






Figure 2: Northern sky as seen by the ARGO-YBJ exper- 
iment in TeV 7-rays. The significance of the excesses in 
standard deviations is shown by the color scale on the right 
side. The two light blue lines indicate the Galactic latitudes 
±5°. 

and AS7, ARGO-YBJ observation on these sources can 
provide more information. Table 2 lists the ARGO-YBJ 
results for the 15 candidate sources, where the Milagro re- 
sult for LAT PSR J2238+59c is given in [19]. 14 out of 15 
sources are observed with significance greater than S.D.. 
9 out of 15 sources are observed with significance greater 
than 1.5 S.D., against 1 expected source from the normal 
Gaussian. The chance probability from Poisson statistics 
would be estimated as 10 -6 . It clearly shows that these 
sources have statistically significant correlations with the 
TeV 7-ray excesses. 

To estimate the energy response and the sensitivity of the 
ARGO-YBJ detector, we simulate different sources in the 
sky with different declinations. Each source is traced by 
means of a complete transit, i.e., 24 hr of observation. The 
simulated events are sampled in the energy range from 10 
GeV to 100 TeV, then the full Monte Carlo simulation of 
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Figure 3: The median energy of 7-ray events that trigger 
ARGO-YBJ and satisfy the event selections as a function of 
source declination. Different lines indicate different spec- 
tral indices from -2.0 to -3.0. 

the detector is performed. Figure 3 shows the median en- 
ergy of 7-rays that trigger ARGO-YBJ and satisfy the event 
selections as a function of the declination of the source for 
several spectral indices. When the index is -2.6, which is 
similar to that of the Crab Nebula, the median energy varies 
from 0.56 TeV at Dec.=30° (the latitude of ARGO-YBJ) to 
2.2 TeV at Dec. =-10° and Dec. =70°, while for sources with 
a hard spectral index -2.0, the corresponding variation of 
median energy is from 1.1 TeV to 4.5 TeV, and the median 
energy varies from 0.36 TeV to 1.1 TeV for sources with a 
soft spectral index -3.0. 

Except the four sources listed in Table 1, the significances 
from all directions are not high enough to claim any sig- 
nal, therefore we set a 90% confidence level (C.L.) upper 
limit to the flux from all directions in the sky. To obtain 
a good data sample with the same detector configuration, 
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Figure 4: 90% C.L. upper limit on the integral flux above 
0.1 TeV averaged on the right ascension direction, as- 
suming a power-law spectrum. Different curves indicate 
sources with different spectral indices -2.0, -2.6 and -3.0. 
The Crab unit is 9.90x 10~ 10 cm~ 2 s^ 1 . 

only events collected with 153 clusters from November 
2007 to February 20 1 1 , are used to estimate the upper limit. 
Firstly, the statistical method given by Helene [20] is used 
to calculate the upper limit on the number of signal events 
at 90% C.L. from each bin. Then the number of events 
is transformed into flux using the full MC simulation tak- 
ing into account the daily path of bins in the sky. Figure 

4 gives the average upper limits to the flux of 7-rays with 
energies above 100 GeV along the right ascension direc- 
tion as a function of declination, which vary between 0.09 
and 0.53 Crab Unit, i.e. 9.90xl0~ 10 cm" 2 s"\ the best 
ones in the world up to now. To investigate the flux upper 
limit for a source with a different spectrum, the upper lim- 
its derived assuming a spectral index -2.0 and -3.0 are also 
presented in Figure 4. The best averaged upper limit for 
sources with index -2.0 is 3% Crab Unit, while the similar 
value for sources with index -3.0 is 20% Crab Unit. The 
upper limits, therefore, depend on the source spectrum. 

It is interesting that Milagro discovered three extended 
sources around the Galactic plane [9], which have flux 
around 1 Crab Unit at 20 TeV. MGRO J 1908+06 and 
MGRO J203 1+41 have been observed by ARGO-YBJ with 
significance greater than 5 S.D., while there is no signif- 
icant excess from the brightest source MGRO J2019+37. 
The upper limit derived using ARGO-YBJ data is lower 
than the flux reported by Milagro even after taking the 
source extension into account, which indicates that the flux 
of this source may not be stable. Details on the analysis can 
be found in [21] . 

5 Conclusions 

The ARGO-YBJ experiment is an air shower array with 
large field of view and has continuously monitored the 
northern sky since November 2007. Using data up to 



February 201 1, we have presented the northern sky survey 
for VHE 7-ray point sources in the declination band be- 
tween -10° and 70°. There are four sources observed with 
significance greater than 5 S.D.. The significance from 
Crab Nebula is about 17 S.D., which indicates that the cu- 
mulative sensitivity of ARGO-YBJ has reached 0.3 Crab 
units. The 90% C.L. upper limits to the 7-ray flux from 
all the directions are obtained under the hypothesis of point 
sources with power-law spectra. The integral flux limits 
above 0.1 TeV vary from 0.09 to 0.53 Crab unit for Crab- 
like source depending on the declination, which is more 
restrictive than previous limits. Meanwhile, ARGO-YBJ is 
still collecting data as the only EAS array at TeV energies. 
The cumulative sensitivity will increase over time. 
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